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We discuss the Josephson current between two noncentrosymmetric superconductors. The coex- 
istence of superconducting order parameters between spin-singlet As and helical p-wave spin-triplet 
At enriches a variety of low-temperature behavior of Josephson current depending on their relative 
amplitudes. We will show that characteristic behaviors of the Josephson current for As > At are 
clearly different from those for As < At. The topologically protected zero-energy surface bound 
states are responsible for the clear difference. We conclude that the Josephson current well re- 
flects character of the topological surface states and the pairing symmetry of noncentrosymmetric 
superconductors. 

PACS numbers: 74.45. -|-c, 74.50. -|-r, 74.25. F-, 74.70.-b 
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I. INTRODUCTION 

Coexistence between the spin-singlet superconducting 
order parameters and spin-triplet one is the essential fea- 
ture of noncentrosymmetric superconductors (NCS)-^— . 
The absence of spatial inversion symmetry leads to spin- 
orbit coupling large enough to mix the spin-singlet com- 
ponent and spin-triplet one. The amplitude of the spin- 
singlet component As and that of the spin-triplet one 
At is a material parameter determined by the amplitude 
of spin-orbit coupling. The Rashba type spin-orbit cou- 
pling induces the helical p-wave spin-triplet order param- 
eter which is the topologically nontrivial superconducting 
state^i^. There have been several studies on supercon- 
ducting properties coexisting of the spin-singlet s-wave 
and the spin-triplet helical p-wave symmetries'^—. It is 
known that topologically protected states with linear dis- 
persion appear at a surface of a NCS for At > As. Re- 
cent papers, however, have suggested a mixed order pa- 
rameter spin-singlet d-wave and spin-triplet p-wave sym- 
metriesi^ii^ which has been proposed for the interfacial 
superconductivity^. Such pairing symmetry results in 
dispersionless surface bound state at the fermi level. A 
similar flat zero-energy surface states has also discussed 
in a NCS very recently^. Physical values originated from 
the bulk region of a superconductor such as the specific 
heat and spin susceptibility'* are expected to be interpo- 
lated from those in the two limits: the pure spin-singlet 
case and the pure spin-triplet one. An open question is 
how physical values governed by the surface bound states 
behave as a function of the relative amplitude between 
At and As- The present paper addresses this issue. 

The surface bound states of unconventional super- 
conductors have been theoretically discussed in heavy 
fermionic superconductors*^, the polar state of "^He*?, 
and high-Tc cuprate o*^i^° . Experimentally, the presence 
of the surface bound states have been observed as the 
zero-bias anomaly^'''^* of the scanning tunneling spec- 
troscopy (STS) of hole-dope d^^i^"^ and electron-doped^^ 
high-Tc cuprates. The zero-bias anomaly has been ob- 



served also in the differential conductance of ramp-edge 
junctions of hole-doped high-Tc cuprates^^ and grain 
boundary junction of electron-doped high-Tc cuprates^. 
The presence of the surface bound states has been rein- 
terpreted since the proposal for new classification of mat- 
ter—. The surface bound states at the zero energy are 
necessary to naturally connect a nontrivial topological in- 
teger number inside of an unconventional superconductor 
with the trivial topological number outside of the super- 
conductor. The dispersion of the subgap states depends 
on the type of the topological number defined in super- 
conductors. The chiral or helical superconductors give 
rise to dispersive surface bound states^^. On the other 
hand, dispersionless zero-energy states are formed under 
dx2_y2- and p^^-wave symmetries. 

In direct current Josephson effect, the surface bound 
states result in large Jci?jv/(Ao/e) values and the de- 
viation of current-phase relationship from the sinusoidal 
function at low temperatur o^^'^" , where Jc is the crit- 
ical Josephson current, Rn is the normal resistance of 
a junction, and Aq is the amplitude of pair potential 
at the zero temperature. Such behavior is called low- 
temperature anomaly of Josephson current and is known 
to be sensitive to spectra of surface bound states^^"'^'^ . So 
far the Josephson effect between s-wave superconductor 
and a NCS has been reported^. The low-temperature 
anomaly of the Josephson current has never been dis- 
cussed yet. 

In this paper, we theoretically calculate the Joseph- 
son current between two NCS's^'^' based on a current 
formula" in terms of the Andreev reflection coefficients 
of junctions. We assume a order parameter which is a 
mixture of the spin-singlet s-wave and the spin-triplet 
helical p-wave symmetries. We call such states as s -I- p 
mixture. For At < As, the Josephson current saturates 
at low temperature as is described by the Ambegaokar- 
Baratoff formula^?. On the other hand for At > As, 
the Josephson current increases logarithmically with de- 
creasing temperature (T). The interfacial bound state 
causes the low-temperature anomaly. The characteristic 



behavior of the Josephson current does not changes grad- 
ually as a function of the relative amplitude between Ax 
and As- The critical point At = As clearly divides the 
qualitative feature of Josephson current. 

In addition to s+p mixture, we also consider two types 
of mixed order parameter between the spin-singlet dxy- 
wave and the spin-triplet helical p-wave symmetries. We 
call such state as d + p mixture. The feature of Joseph- 
son current are well characterized by the relative am- 
plitude between Ax and Ag. In some cases, the Joseph- 
son current increases as 1/T with decreasing temperature 
due to dispersionless zero-energy state. We will discuss 
the physics behind such clear qualitative change of the 
Josephson effect in terms of topologically protected zero- 
energy surface states. It is known that excitation of such 
surface bound states on superconductor are character- 
ized by the Majorana fermior^i^i^. Unusual phenom- 
ena peculiar to the Majorana fermion has been suggested 
theoretically^"—. 

This paper is organized as follows. In Sec. II, we dis- 
cuss a theoretical model of Josephson junction consisting 
two NCS's. In Sec. Ill, we show the calculated results 
of Josephson current for s+p and d + p mixtures. We 
summarize this paper in Sec. IV. 



II. MODEL 

Let us consider a Josephson junction between two 
NCS's as shown in Fig.[l] where the electric current flows 
in the x direction and the junction width in the y direc- 
tion is W. We apply the periodic boundary condition 
in the y direction. The Bogoliubov-de Gennes (BdG) 
Hamiltonian in momentum space reads 
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where (Tj for j = 1 — 3 are the Pauli matrices, ctq is the 
unit matrix in spin space, k^ty) is the wavenumber in 
the x{y) direction, kp is the Fermi wave number, /i is 
the chemical potential, and A is amplitude of the spin- 
orbit interaction. In this paper, we assume that A ^ /i. 
We consider the Rashba type spin-orbit coupling reflect- 
ing the noncentrosymmetry along the z direction (i.e., 
g = {ky,—kx,0)/kp). Correspondingly, we choose the 
d-vector in the pair potential a.s d = g as discussed in 
Ref. 0. As a consequence, the spin-triplet part of the 
pair potential has the helical p-wave symmetry. In this 
paper, we consider three types of mixed order parameter 
as follows 
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where 7 is the incident angle of a quasiparticle as shown 
in Fig. [Hand e'^'' = {k^ + iky)/kp. The first one consists 
of s-wave singlet and helical p-wave triplet components. 
The pair potential of d + p II is the order parameter dis- 
cussed in the interfacial superconductivityi^. Although 
the pair potential oi d + p I may not have a relation 
to any materials, we consider it for theoretical interest. 
The energy eigen values of Eq. ([1]) are E = ±E± with 

E± = JiCk ± A)2 + A^ and 



As±Ax, s+p, 

A±(7)=<( Assin27±Ax, d + pl, 
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To represent the wave function of a quasiparticle, we need 
another values of the pair potential defined by A±{"f) = 
A±(7r-7) 



FIG. 1: (color online). A schematic picture of the Josephson 
junction. 
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The wave function in the left superconductor is obtained as 
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where a± and 6± are the amphtudes of incoming waves, A± and B± are the amphtudes of outgoing waves, and 
Lpj for J = L or i? is the macroscopic phase of a superconductor. In the same way, the wave function in the right 
superconductor is represented by 
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with C± and I?± being amphtudes of outgoing waves. At the junction interface, we introduce the potential barrier 
described by Vq5{x). Throughout this paper, we fix zq = {yam)/{lT?kp) = 5, which leads to the transmission 

probability of the insulating barrier Tb = J^ d'jcoa^ 1 /{zq + cos^ 7) being about 0.01. By eliminating C± and D± 
using a boundary condition, it is possible to obtain the reflection coefficients. 
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The Josephson current can be calculated based on a formula'^^ after applying the continuation E — )■ ia;„ = i(2r7, + l)7rT, 
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a;2 + Aj_ and ri„± 
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A j_ . We introduce a parameter < a < 1 to tune the mixing rate 



between the spin-singlet and spin-triplet components as 

As = aA, At = (1 - a)A, 



(13) 



where the dependence of A on temperature T is described by the BCS theory. The energy spectra of subgap state at 
a surface of NCS is calculated from Eq. (fTO|) with a bounda ry condition '^R{Q,y) = 0. 



III. RESULTS 



At first, we summarize the energy spectra of subgap 
state at a surface of superconductor in Figl2l In s + p 
mixture, energy of bound state satisfies 

(^2 _ A+ A_) cos^ 7 + n+n- (1 + sin^ 7) = 0. (14) 



It has been already known that the surface bound state 
is absent for As > At, whereas the surface bound states 
with the linear dispersion exist for As > At as shown in 
(a). 



In d -|- p I mixture, energy of the surface bound state 
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FIG. 2: (color online). The energy spectra of surface bound 
state, (a): s +p at a = 0, (b): d + p I at a = 0.6, (c): d + p 
I at a = 0.4, (d): d + p II at a = 0.6, and (e): d + p II at 
a = 0.4. The horisontal axis ky/kp corresponds to sin 7. 



satisfies 

£;2(1 + cos^ 7) + sin^ 7(A+A_ + VL+Q.^) = 0. (15) 

The equation has two solutions. The first one is _E = 
which is allowed for |sin7| > Ax/Ag as shown in 
Fig. [2{b). The dispersionless zero-energy bound states 
are a direct consequence of the da;j,-wave symmetry^S. 
Therefore such flat zero-energy state is absent for Ag < 
At. The second solution is given by 



E = ±sin7A/A| - A|4sin^7 



(16) 



for |tan7| < At/(2As) as shown in Fig.[2][b) and (c). 

In d + pU mixture, energy of the surface bound states 
satisfies 

£;2(1 + sin^ 7) + cos^ 7(A+A_ + VL+VL^) = 0. (17) 

The equation has two solutions. The first one is -E == 
for all 7 which is possible only when Ag > Ax as shown 
in Fig. [^fd). The second solution is given by 



E = ±2 cos^ 7a/ A| sin^ 7 - A, 
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which is allowed for | sin 7 1 < -y/Ag/Ax as shown in 
Fig. He). 



In Fig. [3j we show the calculated results of Joseph- 
son current for s -\- p mixture. In Fig. El^a), we plot 
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FIG. 3: (color online). Results for s-l-p mixture. The critical 
Josephson current is plotted as a function of temperature in 
(a) and (c). The current-phase relationships are show in (b) 
and (d) at a low temperature (T = O.OOlTc). The transmis- 
sion probability of the tunnel junction Tb is fixed about 0.01 
throughout this paper. 

the critical Josephson current as a function of temper- 
ature for several choices of a satisfying Ag > At. The 
Josephson current is normalized by Jq = 7rAo/(2ei?Ar), 
where Rn is the normal resistance of the junction. In the 
case of metallic superconductor junctions, the Joseph- 
son critical current becomes Jo at the zero tempera- 
ture. In Fig[3jb), we show the current-phase relation- 
ship (CPR) at a low temperature T — O.OOlTc for a = 1, 
0.8 and 0.6, where (f = ipL — fR is the phase differ- 
ence across the junction. The critical Josephson cur- 
rent saturates at low temperature and the CPR is si- 
nusoidal for a — 1, 0.8 and 0.6. Thus the Josephson 
current for Ag > At obeys the Ambegaokar-Baratoff 
relation because there is no surface zero-energy states. 
In Fig. ^c) and (d), we respectively show the depen- 
dence of critical current on temperature and the CPR 
at a low temperature for several choices of a satisfy- 
ing Ag < At. The results in Fig. [Sfc) and (d) show 
qualitatively different behavior from those in Fig. El^a) 
and (b), respectively. The critical Josephson current 
in Fig. [3fc) increases with decreasing temperature even 
far below Tc- This behavior is called low-temperature 
anomaly of Josephson current. The resonant tunneling 
through the surface bound state at the zero energy is re- 
sponsible for the anomal y ^^1'^° . Such zero-energy state is 



possible sin 7 = as shown in Fig. [2][a) . According to 
previous paper a^^'^^ , the Josephson critical current in- 
creases logarithmically with decreasing temperature for 
a = 0. The results for a = 0.2 and 0.4 also show the 
logarithmic low-temperature anomaly. Correspondingly, 
the contribution of higher harmonics slightly deviates the 
CPR from the sinusoidal relation as shown in Fig. [2Id). 
Thus the characteristic feature of Josephson current qual- 
itatively changes at the singular point of As = At. 




B. d+p I 

Next, we show the calculated results of Josephson cur- 
rent for d + p 1 mixture as shown in Fig. U The tem- 
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FIG. 5: (color online). The current-phase relationship in Fig. 
4(b) and (d) are calculated for a higher temperature at T = 

o.ir,. 



temperature anomaly. The presence of the zero-energy 
surface bound state at sin 7 = in Fig. [SJc) explains 
the similarity. Thus the characteristic feature of Joseph- 
son current for d + pl mixture also qualitatively changes 
around the point of As = At- 

The large deviation of CPR from the sinusoidal func- 
tion in Fig. m^b) can be seen only at low temperature. 
According to an analytical expression of Josephson cur- 
rent for a = 1, the higher harmonics contribute to the 
Josephson current when the temperature is much smaller 
than v^TbAq. Here Tb is the transmission probability 
of the tunnel junction and is about 0.01 in the present 
calculation. In Fig. [SJ we show CPR for a higher tem- 
perature at T = O.lTc. In both Figs. Hlja) and (b), the 
CPR deviates from the sinusoidal relation only slightly 
at T = O.lTc. On the other hand, the amplitudes of the 
Josephson current remain sufficiently larger value than 
Jo. 



C. d+p II 



FIG^ 4: (color online). Result for d + pl. The critical Joseph- 
son current is plotted as a function of temperature in (a) and 
(c). The current-phase relationship is show in (b) and (d) at 
low temperature T = O.OOlTc. 

perature dependence of critical current for a =1, 0.8, 
and 0.6 satisfying As > At show the low-temperature 
anomaly as shown in Fig. Sl^a). The critical current in- 
crease as T~^ with decreasing temperature at a = 1—. 
The results for a = 0.8 and 0.6 also show such power- 
law like low-temperature anomaly. Corresponding CPR 
shown in Fig. Sljb) indicates jump at 1^9 = tt because 
of the contributions of higher harmonics. In this case, 
the surface bound state are energetically localized at 
■B = as shown in Fig. [Hb). The resonant tunnel- 
ing through such zero-energy states causes the strong 
low-temperature anomaly. For As < At, on the other 
hand, the Josephson current in Fig. Hl^c) and (d) have 
properties similar to those in Fig. ^c) and (d), respec- 
tively. The critical current indicates the logarithmic low- 



Finally, we show the calculated results of Josephson 
current for d + p 11 mixture as shown in Fig. [51 The 
temperature dependence of critical current for several a 
satisfying As > At indicate the strong low-temperature 
anomaly as shown in Fig.[6lja) and (b). The critical cur- 
rent increase with decreasing temperature as T~^ and 
the CPR at a low temperature shows the jump at ip — n. 
The zero-energy surface bound states are possible for 
all 7 as shown in Fig. [Ifd). The presence of the flat 
zero-energy states explains the similarity of the results 
in Fig.inia) and (b) to those shown in Fig.|3l^a) and (b), 
respectively. The calculated results for a =0.4, 0.2 and 
satisfying As < At are shown in Fig. IBJc) and (d). 
On the contrary, the results for As < At in Fig. ^c) 
and (d) has properties similar to those in Fig. [2Ia) and 
(b). Namely, the Josephson current saturates at low tem- 
perature and the CPR is sinusoidal at low temperature. 
The zero-energy state at sin 7 = ±1 exists as shown in 
Fig- Efe). Although this zero-energy state appears as a 
results of the resonant Andreev reflection^^, it does not 



so much affect the Josephson current. The wavenum- 
ber sin 7 = ±1 means k^ = cos 7 = 0. Thus a quasi- 
particle does not have momenta in the current direction 
in the zero-energy state. When we consider huge spin- 
orbit coupling, it has been pointed out that^^i^i^ the 
flat zero-energy states appear for sin7 > (1 — 2A//i). 
In such case, the flat zero-energy state may cause the 
low-temperature anomaly. This statement, however, is 
still unclear in realistic junctions with a thick insulat- 
ing barrier because the contribution of a quasiparticle 
with sin 7 « ±1 to Josephson current becomes expo- 
nentially small. Within the approximation of A//i ^ I, 
there is no effective zero-energy state which causes the 
low-temperature anomaly for Ag < At. Therefore the 
Josephson current in Fig. |6l[c) and (d) show qualitatively 
the same behavior as those in Fig. [31[a) and (b), respec- 
tively. Thus the characteristic feature of Josephson cur- 
rent for d + p 11 also qualitatively changes around the 
point of Ag ~ At. 

At Ag = 0, subgap state with the linear dispersion 
appears around sin 7 = as mathematically shown in 
Eq. p^ . This zero-energy state, however, is not a result 
of the resonant Andreev reflectionSi but is a result of a 
node in the pair potential. Thus Ag = cannot be a 
critical point. Actually, the results for a = in FiglHl^c) 
show the saturation of the Josephson critical current at 
low temperature. 
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FIG. 6: (color online). Results for d-l-p II. The critical Joseph- 
son current is plotted as a function of temperature in (a) and 
(c). The current-phase relationship is show at low tempera- 
ture r = 0.001Tc. 



IV. CONCLUSION 



In summary, we have theoretically studied the Joseph- 
son current between two noncentrosymmetric supercon- 
ductors based on the Bogoliubov-de Gennes equation 
and a general current formula. We have assumed three 
types of order parameter which consists of spin-singlet 
Ag and spin-triplet components At at the same time. 
The Josephson current for Ag > At shows clearly the 
different characteristic behavior from those for Ag < At 
for all pairing symmetries. The clear difference can be 
understood by analyzing the topologically protected zero- 
energy states at a surface of noncentrosymmetric super- 
conductor. The dispersionless zero-energy bound states 
are responsible for strong low-temperature anomaly of 
Josephson current in which the Josephson critical cur- 
rent increases as l/T with decreasing temperature. The 
surface state with linear dispersion causes the weak low- 
temperature anomaly in which the Josephson critical cur- 
rent increases logarithmically with decreasing tempera- 
ture. When the surface zero-energy state is absent, the 
Josephson current obeys the Ambegaokar-Baratoff for- 
mula. 
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